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Abstract 

We calculate constraints on radiatively decaying neutrinos from the recent detection of singly ionized helium 
in the diffuse intergalactic medium (IGM) at z ~ 3.3. We consider a model in which neutrinos predominantly 
decay into invisible relativistic particles with a rate and with a small branching ratio into the radiative 
mode. To satisfy the observation of singly ionized helium, which puts a lower bound on the number density of 
singly ionized helium in the IGM, we show that: for r ^ 10^® sec(l eV/m^), transition moment of neutrinos 
/ri 2 is constrained to be < 4-8 x 10 “^”^/Tb for 110 eV ^ TOj. ^ lOkeV. We compare this bound with other 
astrophysical and cosmological bounds on radiatively decaying neutrinos. 

PACS numbers: 98.60.Hj, 98.80.Gq, 98.70.Vc, 98.80.bq 


Radiatively decaying neutrinos have found widespread application in cosmology and astrophysics [1- 
3]. On the other hand, a formidable array of constraints exist on radiative lifetime of radiatively decaying 
neutrinos, from studies of various disparate phenomena in cosmology and astrophysics [3-8]. In this paper, 
we point out another constraint which comes from the existence of primordial elements (hydrogen and 
helium) in the diffuse IGM at high redshifts. The search of primordial elements in the diffuse IGM is crucial 
to understanding the universe at early epochs (Gunn-Peterson (GP) test [9]). So far there is no definite 
evidence of the presence of either neutral hydrogen or neutral helium in diffuse IGM [10,11]. Last year, first 
detection of singly ionized helium in diffuse IGM was reported by Jakobsen et al.. [12]. Jakobsen et al. 
detected strong absorption from singly ionized helium along the line of sight of quasar Q0302-003 at z ~ 3.3. 
Their 90% upper bound on the continuum optical depth from singly ionized helium was > 1.7. As there 
are only upper bounds on r^p and T^p\ the presence of singly ionized helium in the IGM puts stringent 
constraints on the magnitude and spectrum of homogeneous sources of photoionization [13]. Radiatively 
decaying neutrinos is one such homogenous source of photoionization. In this scenario, a massive neutrino 
(lyfi) decays into a photon ( 7 ) and a light neutrino (t'^) [14,15]. If m^j, <C —the case we consider—the 
photon and the light neutrino equally share the rest energy of i^h ■ We assume the heavy neutrino to be tau 
neutrino while the light neutrino can be either or ly^. To ionize singly ionized helium the mass of tau 
neutrino must be greater than 108.8 eV (the ionization potential of singly ionized helium is 54.4 eV). The 
closure density constraint restricts the mass of a stable neutrino (or a neutrino which decays with a rate less 
than to being the present age of the universe) to be less than 91h^ eV. This constraint forces the rate 
of neutrino decay to be greater than tg^. As a radiative decay lifetime of less than tg is forbidden by several 
constraints, the neutrinos must predominantly decay into particles which are invisible (majoron, neutrinos 
etc.) with decay rate ^ tg^, with a small branching ratio B for the radiative mode (radiative lifetime 
Tj = ■ 

In the presence of decaying neutrino, the thermal history of the universe can change considerably. First, 
as soon as the neutrinos become nonrelativistic the universe becomes matter dominated (z = Zi). The 
universe remain matter dominated, with massive neutrinos dominating the energy density, until the time 
the neutrinos decay (z = Zd) into relativistic decay products. (For simplicity we assume that all relativistic 
products of decaying neutrinos remain relativistic upto the present epoch.) Thus the epoch of neutrino 
decay is followed by a period of radiation domination in the universe. Depending on the mass and lifetime 
of decaying neutrinos and the energy density of other nonrelativistic matter (Djvr) in the universe, another 
period of matter domination could follow (z = 22 )- In the foregoing we have assumed the neutrinos to 
decay instantanously, which would result in a sudden change in matter and radiation domination. It is a 
good assumption for calculating the time-redshift relationship and for computing the age of the universe 
in the presence of decaying neutrinos [16]. In this paper we assume the universe to be spatially flat i.e., 
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r^o — + ^NR = Ij where ^Inr is the energy density due to all other non-relativistic matter. 

Helium is the second most dominant primordial element ( 8 % by number). In a uniform IGM, the 
proper density of helium is 

UHe = 6.8 X cm“^. (1) 


Throughout this paper we take fl/cM = 0.05 and h = 1/2. For a uniform IGM, any resonance line at 
wavelength A and oscillator strength f, of singly ionized helium with a proper number density nueuiz), will 
produce an optical depth [17,18] 

rS;" = —fH-\z)Xnn.u{z). (2) 

TTleC 

The dominant scattering is at wavelength 304 A. Using this, > 1.7 at z = 3.3 translates into a lower 
bound on the fraction of singly ionized helium, y 2 > 1.2 x 10“^, where j /2 = nHeii/nHe- Using that lower limit 
on j /2 is much less than unity, equilibrium between recombination and ionization processes can be assumed, 
which allows one to express 7/2 as 

7/2 0.026 X 10-1°(3) 

ny{z) 


Here 


mu/2 

/ -jrfdik,t)exp{-T{z,Ze,k) cm 
€0 
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Co = 54.4 eV is the ionization potential of singly ionized helium; we have assumed the temperature of IGM 
T = 1.5 X 10^ K; fd{k, t) is the distribution function of decay photons; t(z, Ze, k) is the optical depth (to be 
described below) suffered by a photon which is emitted at Ze and observed at z with a wave number k. In 
the limit of small energy transfer between electrons and decay photons (a valid assumption for Zd ^ lO"!), 
the distribution function of photons can be written as [1,8] 


fdik,t) = 


BTr'^nvir) 


n{miyj2y^'^k^ it 


3n-l 


exp 


— (t/r) (2A:/mi/)^^" Q{m^/2 — k). 
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Here n = 1/2, 2/3 for radiation-dominated and matter-dominated epochs respectively. 

Before discussing the results of our analysis, we briefly review some of the other bounds on radiatively 
decaying neutrinos. If neutrinos have a coupling to photons then, during a supernovae explosion, a photon 
flux coincidental with the neutrino flux from the supernovae should be observed. No such flux was observed 
during supernovae SN1987A. This has been used to constrain the parameter space of radiatively decaying 
neutrino [6]: 

TiTii, > 2 X 10i°i? eV sec, (6) 


for neutrino masses between 100 eV and a few MeV. Other major astrophysical constraint comes from 
the stndy of red giant stars. The dominant energy loss mechanism in red giants is the decay of plasmons, 
formed in the core of the star. The rate of plasmon decay is enhanced if neutrinos have a magnetic moment. 
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This can lead to precipitous cooling of the star, which is in conflict with known stellar evolution time scales. 
A study of such a process constrains the neutrino transition moment /ii 2 ^ which translates to [7] 


r ^ 2.1 X lO^^Bsec 


/ \ 

VleV/ 
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for rrii, < 10 keV. 


( 6 ) 


An interesting bound on trasition moment of radiatively decaying Dirac neutrinos can be obtained if the 
universe had a primordial magnetic field. The presence of magnetic field can cause the right-handed neutrino 
to be in theraml equilibrium with other relativistic particles in the universe at the time of nucleosynthesis, 
thereby adding an extra relativistic species, which is unaccetable. From such an analysis, a lower bound on 
trasition moment /ii 2 < 5 x 10“^"‘(30/rs can be obtained [19]; or in terms of neutrino lifetime 


T > 2.5 X 10^2 


/30eV\ 

\ rn„ J 


4/3 


sec. 


However, as there is no evidence of the existence of primordial magnetic field, this constraint may not apply. 

Another stringent constraint on radiatively decaying neutrinos comes from the planckian nature of CBR. 
Photons emitted by decaying neutrinos can heat the electrons in the IGM; hot electrons can transfer this 
energy to CBR photons. If this process takes place at redshifts ^ 10®, the spectrum of CBR is altered [20]. 
As CBR is known to be a planckian to a very high accuracy [21], the decay of neutrinos into photon mode 
is severly restricted [8] 

Finally, the decay photons constitute a diffuse background at wavelengths ^ rriyl^. In recent years, 
rocket-borne experiments have placed upper bounds on diffuse, background photon flux in ultra-violet (UV), 
a wavelength range in which radiatively decaying neutrinos contribute most significantly [22]; in addition 
there exist several bounds on the extragalctic hydrogen-ionizing flux from astrophysical observations [23]. 

In this paper we point out that radiatively decaying neutrinos with masses above 108.8 eV are subject 
to an additional constraint: the existence of singly ionized helium in the ICM. From eq.(3), it is clear that if 
n^{3.3) ^ 10“”^, observed bound on r"*” would be violated. To calculate n~^{z), sources of absorption must 
be identified. We consider two sources of absorption: diffuse ICM and lyman-a systems. Optical depth due 
to diffuse ICM is 

UHer'’HeIl(fc)j/2) (8) 

where CTHeii(fc) is the photoionization cross section from the ground state [24]. Tdoud, the avarage attenuation 
of photon flux due to poisson-distributed clouds, is given by 

pz roo 

Tc\ond{z,Ze,k) ~ / / dz d/Vni T’(A^Hi, 2:){ 1 “ exp[-/VHeiiCTHeii(A:)]}. (9) 

J Ze -/O 


Tdi!!nse{z,Ze,k) CZ / dZ 
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Here A^hi and -/Vneii correspond to column densities of neutral hydrogen and singly ionized helium in the 
clouds and V{Nia^z) is the number density of clouds in a given column density and redshift interval. In 
ionization equilibrium, IVneii can be inferred from known column densities of neutral hydrogen A^hi 

n™ 

^He„=iVmX^ Cm-2; (10) 

here n™ is the effective number of photons which ionize neutral hydrogen. From ’proximity effect’ [25], n™ 
is known to approximately constant for 2 ^ 2 ^ 4 with a lower limit n™ ^ 6.2 x 10“®cm“^. For some of 
the paramter space for which ^ 1, the value of n™ from decaying neutrinos is less than the lower limit 
from ’proximity effect’. Hence we include the contribution from ’proximity effect’ in n™ . It is necessary to 
include this additional source of photons to correctly determine the optical depth (eq.(lO)). For V{N, z), we 
take models A1 and A2 of Miralda-Escude and Ostriker [18], in the redshift range 1 < z < 5. 

Our results are shown in figs. 1-4, along with other constraints. To gauge the effect of absorption, 
results in case of no absorption are also shown. In some range of masses, constraints on B from the detection 
of singly ionized helium can be a few orders more stringent than any other constraint. Our results can be 
qualitatively understood. The photon flux at any energy E is determined by the the shape of decay photon 
spectrum (eq(5)). At any given redshift, if the exponent in photon spectrum can be neglected for photon 
energies ^ 54.4 eV then photon flux at these energies 

T{E) (xE^ X fd{E) ~ 


Thus photon spectrum is extremely hard (for comparsion, J- cx E~°‘, 0.7 < a < 1.5, for a photon background 
dominated by quasars). The bounds on B are most stringent in this region. As rup is decreased for a given 
T, the exponent in decay photon spectrum begins to be important. This leads to a decrease in the number of 
photons above energy 54.4 eV as well as to an increase in absorption (eq(lO)), because the photon spectrum 
becomes soft for energies ~ 54.4 eV. And so larger values of B are required to ionize the singly ionized 
helium. Our results are valid for 10^^ sec < t < 10^® sec. It is because for r much less than 10^^, the 
assumption that interaction between photons and electrons can be neglected in determining the spectrum of 
the decay photon breaks down. For r > 10^® sec, neutrino masses above 108.8 eV violate the closure density 
constraints. (All the region shown in figs (l)-(4) is compatible with closure density bounds.) The bounds on 
radiative lifetime can be translated into bounds on magnetic moment of radiatively decaying neutrinos. 
For Dirac neutrinos, radiative lifetime Tj (= t/ B) is related to neutrino transition moment /ii 2 as [26] 


= 10 ^® 


io-^Vb 


The bounds on /ii 2 from the existence of singly ionized helium in the IGM can be summarized, for 
model A2 for absorption from lyman-a systems, as: 
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(1) For r > 1.6 X 10^® sec(l eV/m^), 

Mi2 ^ 4-8 X 10“^”^/is, 

in the mass range llOeV ^ ^ lOkeV. 

(2) For 1.6 X 10^®sec(leV/TO^) ^ r ^ 5 x 10^^ sec(l eV/m^), 

^il2 ^ 8-100 X 10“^”^ /is, 

in the mass range llOeV ^ ^ lOkeV. 

For T ^ 5 X 10^"^ sec(l eV/m^), other bounds become more signihcant (Hgs. 1-4). 

Our bounds on /ii 2 are not restrictive enough to constrain the standard Wienberg-Salam model. How¬ 
ever, several extensions of the standard model like the left-right symmetric model, charged-Higgs model, and 
broken R-parity supersymmetric model can give large transition moments (for a recent review see [27], [28]), 
and therefore would be constrained by the bounds we derive. As the most stringent constraints are obtained 
in the range r ^ 1.6 x 10^® sec(l eV/m^), our results can be used to put an upper bound on t. It should 
be pointed out that the though these values of r are not constrained by the closure density constraints, age 
and structure formation constraints may rule out much of the range of r in which most stringent constraints 
on fj,i 2 are obtained ([29], [30]). 

We now discuss various sources of errors in our estimates from uncertainties in input parameters. (As 
GP tests are quite sensitive and model independent, our analysis is free from any errors of modelling.) First, 
it has been pointed out that a large contribution to Jakobsen et al.’s bound of r“p” > 1.7 might come from 
line blankating [13]. But the presence of a foreground quasar along the line of sight to Q302-003 might 
reduce the optical depth due to line blankating quite significantly [31]. In any case, these uncertainties 
should not change bounds on B by more than an order. (It should however be pointed out that even in 
case most of the contribution to optical depth observed by Jakobsen et al. comes from line blankating, 
equally stringent bounds on radiatively decaying neutrinos can be obtained. If line blankating causes most 
of the optical depth, one requires Sl = ^ 65. This puts strong constraints on the spectrum of 

radiatively decaying neutrinos.) Other input parameters like flioM, T, and h can change our estimate by 
a factor of few. Another source of uncertainty lies in the choice of a model for absorption from lyman-a 
system. From observations, Model A2 of Miralda-Escude and Ostriker gives a fair picture of absorption due 
to lyman-a systems at high redshifts. However, Meiksin & Madau [32] have argued that absorption due to 
lyman-a clouds can be smaller (as compared to model A2 of Miralda-Escude and Ostriker). On the other 
hand, Lanzetta [33] has reported a rapid evolution of lyman-limit systems at high redshifts, which would 
increase the absorption. However these uncertainties in absorption from lyman-a systems do not change 
values calculated using model A2 by more than a factor of few. 

To conclude, we showed that the presence of singly ionized helium in the diffuse IGM can put extremely 
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stringent constraints on radiatively decaying neutrinos. In particular, for r ^ a few x 10^® sec(leV/m^), the 
bounds on radiative lifetime (or equivalently on transtion moment ^ 12 ) can be several orders more stringent 
than the previously known astrophysical and cosmological constraints. More recently, large optical depth 
from singly ionized helium has been detected in the spectra of two more quasars [34], which would allow one 
to obtain similar constraints on radiatively decaying neutrinos. Our results also show that the presence of 
primordial elements in the IGM generically give very restrictive bounds on radiatively decaying neutrinos. 
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Figure Captions 


Figure 1 


Figure 2 
Figure 3 
Figure 4 


: For T = 10^^ sec, constraints from singly ionized helium (solid lines marked with chosen model of 
absorption from lyman-o; systems) are plotted with constraints from supernova (dashed line), CBR 
spectrum (dotted line), UV background for absorption model A2 (dot-dot-dot-dashed line), and primor¬ 
dial magentic field (dot-dashed line). Region above the curves is ruled out by the respective constraints. 

: Same as fig. (1) but for t = 10^^ sec. 

: Same as fig. (1) but for r = 10^^ sec. 

: Same as fig. (1) but for t = 10^® sec. 
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